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ABSTRACT 

Nuclear reactions In the hot accretion plasma surrounding a collapsed 
star are a source of neutrons, primarily through spallation and pion-producing 
reactions, and antineutrons, principally through the reaction p+p+p+p+n+n. We 
calculate spectra of neutrons and anti neutrons produced by a variety of 
nonthermal energetic particle distributions in which the target particles are 
either at rest or in motion. If only neutral particles are free to escape the 
interaction site, a component of the proton and antiproton fluxes in the 
cosmic radiation results from the neutrons and antineutrons which leave the 
accretion plasma and subsequently decay in the interstellar medium. This 
additional antiproton component could account for the enhanced flux of 
antiprotons in the cosmic radiation, compared to values expected from the 
standard "leaky-box" model of cosmic-ray propagation and confinement. 

Moreover, the low-energy antiproton flux measured by Buffington et al. (1981) 
could result from target-particle motion in the accretion plasma. This model 
for the origin of antiprotons predicts a narrow 2.223 MeV line which could be 
observable. 
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INTRODUCTION 

Matter accreting onto collapsed objects can attain very high energies due 
to the conversion of gravitational energy into kinetic energy. Because of 
thermal decoupling of the electrons and ions in accreting thermal plasma, the 
ion component may reach temperatures In excess of 100 MeV in disk (e.g., 
Shapiro, Lightman and Eardley 1976; Eilek and Kafatos 1983) and spherical 
accretion models (e.g., Dahlbacka, Chapline and Weaver 1974; Kolykhalov and 
Sunyaev 1979; Meszaros and Ostriker 1983). Neutrons are formed in a hot 
accretion plasma through nuclear breakdown reactions at temperatures > 1 MeV 
(Aharonian and Sunyaev 1984; Gould 1986) and In association with pion- 
producing reactions (Dermer 1986) at temperatures > 20 MeV. At temperatures 
> 200 MeV, a channel of neutron and anti neutron production is opened through 
the reaction p+p+p+p+n+n which, moreover, produces a spectrum of 
antineutrons extending to low energies (Dermer and Ramaty 1986), in contrast 
to the spectrum of antineutrons produced In collisions of protons with a 
stationary target (e.g., Stephens 1981). These neutrons and antineutrons, 
unconfined by the magnetic field of the accretion plasma, can escape into the 
interstellar medium and, after decaying, be detected as a low energy proton 
and antiproton component of the cosmic radiation. 

Support for thermal accretion models of galactic X-ray sources is 
provided particularly by the case of Cyg X-l, whose X-ray spectrum is 
successfully fit using either disk (Shapiro, Lightman and Eardley 1976) or 
spherical (Meszaros 1983) accretion geometries. But the existence of galactic 
point sources emitting gamma radiation at energies e>10 12 eV implies that 
energetic nonthermal particle spectra are common in galactic binary X-ray 
sources such as Cyg X-3 and Her X-l. Particle acceleration in confined 
volumes, with maximum sizes inferred from time-variability arguments, requires 
intense magnetic fields (Hi 11 as 1984a). Only charged particles with energy 
E(GeV) > 90ZB 3 M o , where B 3 is the magnetic field at the Interaction site in 
units of 10 3 gauss, M Q is the mass of the collapsed object in solar masses and 
Ze is the charge of the ion, can directly escape from the interaction region, 
assuming that the gyroradius must approximately equal the Schwarzschild radius 
for charged-particle escape. Neutrons and anti neutrons, on the other hand, 
leave the interaction region directly provided that (i) the mass of the 
collapsed object is <10 M 0 (Dermer and Ramaty 1986), (ii) the grammage 
shrouding the interaction region is < 50 gm/cm 2 , and (iii) u*3B/dx » 
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10 _23 B 3 /M o < E(GeV), where y is the magnetic dipole moment of the neutron. 
These conditions ensure that the neutrons or anti neutrons (i) do not decay 
before escape, (11) do not Interact strongly before escape, and (ill) do not 
significantly change energy when traversing the intense, turbulent magnetic 
field likely to be found in the accreting plasma. Condition (iii) is easily 
satisfied even if the correct distance scale for magnetic field variation is 
much smaller than the Schwarzschild radius used in this estimate. The other 
conditions are also not difficult to satisfy in galactic sources. 

In this paper, we calculate the spectra of neutrons and antineutrons 
produced by energetic nonthermal particle distributions located in the 
vicinity of a compact object. We do not consider the detailed nature of the 
acceleration mechanism and interaction site, but for simplicity assume an 
Isotropic, energetic particle distribution of solar (photospheric) 
composition. We consider two interaction models which take into account the 
trapping of the charged particles in the interaction region. In the first, 
the particles are assumed to Interact with an ambient background plasma in the 
thick-target model (e. g., Ramaty, Kozlovsky and Lingenfelter 1975). This 
model has been extensively developed to Interpret gamma-ray and neutron 
observations from solar flares (Murphy and Ramaty 1985, and references 
therein). In the second, the energetic particles interact only with each 
other, that is, we consider a relativistic nonthermal plasma. We employ 
shock-acceleration spectra with a spectral index -2.2 appropriate to the 
cosmic radiation before effects of the energy-dependent escape are taken into 
account. This source could account for the antiprotons observed in the cosmic 
radiation (Golden et al. 1979; Bogomolov et al. 1979; Buffington et al. 1981) 
in excess of the leaky-box model prediction. A narrow 2.223 MeV line due to 
the energetic neutrons thermal i zing and being captured by protons in the 
atmosphere of a binary stellar companion might also be formed. 

ANALYSIS 

In the thick-target model, the interacting primary and secondary 
energetic charged particles remain trapped in the interaction region and lose 
energy through particle collisions with the ambient plasma, while the neutral 
particles are free to escape. The spectrum of the antineutrons injected into 
the interstellar medium is identical to the spectrum at production under the 
conditions listed in the Introduction, except for a possible gravitational 
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redshiftlng. Neglecting this last complication, the source spectrum of 
antineutrons (particle 3) produced in the reaction 1 + 2 ♦ 3 + X is just 


n 3 (E 3 ) [cm" 3 -s~ 1 -GeV _1 ] 


4irnJdE 2 .J 2 (E 2 ) 

E? 1n (E,) 


d °l+2-3+X^ E 3i E 2) 


dE, 


( 1 ) 


where nj is the ambient target particle density, J 2 (E 2 ) is the equilibrium 
flux of energetic particle species 2, including secondaries, E 2 1n (E 3 ) is the 
minimum energy of a particle of type 2 necessary to produce particle 3 with 
energy E 3 , and do^ +2+3+x (E 3 ;E 2 )/dE 2 is the differential cross section for the 
production of particle 3 with energy E 3 in a collision between particle 2, 
with energy E 2 , and stationary particle 1. The flux [cm _2 -s' 1 -sr* 1 -GeV _1 ] 

E‘ 

dE ® 2 dE 

0j(E 2 ) - [4*11,1— I r l J’dEJ-n 2 (EJ)-exp(-J'dE“/[A 12 (E“)- 1 — -| 1 } , (2) 

H dA E 2 E 2 dA 

where n 2 (E 2 ) is the rate at which particles of type 2 are injected per cm 3 
per GeV, and dE 2 (E 2 )/dA is the rate at which they lose energy per unit 
grammage traversed. Nuclear collisions are treated as a particle sink in 
terms of the grammage a 12 (E 2 ) » mp/oi 2 (-E 2 ), in which case secondaries 
particles must be added explicitly in the term n 2 (E 2 ) (see below). 

In a relativistic plasma, the energetic particle density is sufficiently 
great that reactions of particles in motion with other particles in motion 
dominates secondary production. A relativistically correct description of the 
reaction rate for producing particles of type 3 in reactions of particles of 
type 1 and 2 is given by 

3 1 l “ 00 + 1 

Q 3 lcm" -s' ] = %(1+« 12 )“ cf dp iF 1 (p 1 )/ dp 2 F 2 (p 2 )Jd cose-s • (l-e^cose) 

0 o _1 ' 

x fe’ 3 d2 °Wx( E *-^=v) (3) 

m 3 -1 dE*dy* 


(cf. Dermer 1984). In this equation. 



(4) 


is the density, F-j(P-j) is the isotropic differential momentum spectrum, and 
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Is the momentum variable of particles of type i. The term 
E* max 3 1S the maximum CMS energy of particle 3 and is a function of the 
relative Lorentz factor Y r =(E 1 E 2 -p 1 P 2 cose)/m 1 m 2 onl y» where E x and E 2 are the 
energies of the colliding particles and 0 is the angle between their 
directions of motion in the LS. The CMS differential cross section for the 
production of particle 3 in terms of CMS energy E* and CMS angle e*=cos _1 p* is 
denoted by d 2 a*/dE*dp*. 

We employ nonthermal energetic particle spectra predicted by shock 
acceleration theory. The momentum spectrum resulting from diffusive shock 
acceleration is given by F^p^) « p^ (e.g., Blandford and Ostriker 1978). 
Normalizing this spectrum using equation (4) gives 


n i = dp r pT q , 1-1,2, 

p co,i 


(5) 


where Kj is the normalization constant, n^ is the volume density, and p COjl - is 
the low-momentum cutoff in the momentum spectrum of particles of type i. 


We solve equation (3) using the Monte-Carlo method, generalizing the 
technique of Ramaty and Meszaros (1981) to secondary particles created in the 
collisions of particles with, in general, different masses. Introduce an 
array of six numbers chosen randomly between 0 and 1, namely r!,...,r 6 . 

Choose values of momenta p^ and P2 from the inverted particle distribution 
functions. The shock spectrum can be analytically inverted, giving 


Pi = P, 


<Fq> 


CO,l 1 


1 - 1 , 2 . 


( 6 ) 


Choose the remaining variables from uniform distributions: cose=2r 3 -l, 

E* = m 3 +[ E* ax 3 (Y r )-m 3 ] *r 4 , p*=2r s -l and <t»*=2itr 6 , where <t>* is the CMS 
azimuthal angle of the outgoing particle of mass m 3 . The reaction rate (3) is 
then given by 


N 


2S y 
r r i 


%n.n 2 c "tot 
q 3 - » tot (i+« 12 > ,1 


dE*du* 


(7) 


in the limit N tot , the number of Monte-Carlo trials, goes to infinity. 
Differential energy spectra are determined by calculating differential 
reaction probabilities from the individual Monte-Carlo trials of equation (7) 
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and binning in energy. The energy of particle 3, E 3 , was determined from the 
kinematical relations E 3 = Y c (E* + 0 c P*cose*), where y c = (yx+YzJ/s* 5 * 

2 2 2 %c 2 

s = m 1 +m 2 +2m 1 m 2 Y r , and cos e* = y*n*+(l-u* ) -(1-n* ) -cos <t>*, 

2 2 J* 

where n* - ( e i-y c E*)/s c y c P* and Ef = (s+m 1 -m 2 )/2s^. In the limit Y 2 -*-l» that 

is, F(p 2 ) -► 6 (p 2 ) , equations (3) and (7) reduce to the stationary target 
case. 


NUMERICAL RESULTS 

A detailed description of the cross sections and energy spectra for the 
production of secondary neutrons and protons in reactions involving protons 
and a particles is given by Murphy, Dermer and Ramaty (1987). A description 
of secondary neutron production in collisions involving heavier nuclei will be 
given elsewhere. A complete version of the calculations relevant to this work 
is in preparation. 

For calculations of neutron and antineutron production in a thick target, 
we assume solar composition (Cameron 1982) for the ambient plasma and the 
energetic particles, that is, the number of particles with energy greater than 
30 MeV/nucleon is assumed proportional to the Cameron abundance of that 
species. Calculations of secondary proton and neutron production at energies 
< 5 GeV/nucleon involve integrations over complicated secondary energy spectra 
and cross sections and were evaluated from equations (1) and (2) using a Monte 
Carlo simulation (cf. Murphy 1985). At energies > 5 GeV, the effects of a 
particles and heavy nuclei can be treated as a multiplicative correction to 
the basic proton-proton (p-p) reactions. Taking A p _ p = 50 gm/cm and | dE p /dA | 
= 4xl0 -3 GeV/(gm/cm 2 ) (Murphy 1985) in equation (2) gives a simple 1- 
dimensional integral for the thick-target equilibrium proton flux. Using this 
flux, we calculated the injection sources of anti neutrons and > 5 GeV neutrons 
from equation (1) using the scaling representations of Tan and Ng (1983a; see 
Murphy, Dermer and Ramaty 1987). 

In the calculations of neutron and anti neutron production in a 
relativistic plasma, we assume that there is 1 proton per cm 3 with momentum 
greater than the low-momentum cutoff, and that the number of energetic 
particles above the cutoff momentum is proportional to the Cameron abundance 
of that species. If the low-momentum cutoff p CQ ^ tends to zero, interactions 
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of the energetic particles with particles nearly at rest dominate secondary 
production because of the divergence in the normalization of the shock 
spectrum in this limit [see eq. (5)]. The problem therefore reduces to a 
thin- target calculation. As P COt -j becomes large, the average collision energy 
increases, and reactions of particles in motion with other particles in motion 
becomes more probable. Due to the large antineutron production threshold, the 
effect of increasing P co j causes the relative production of antineutrons 
compared to neutrons to increase. Also, the spectrum of anti neutrons now 
extends to low energy, as most n are formed with low energy in the CMS of the 
colliding protons, and the CMS can now have low velocity in the rest frame of 
the plasma (Dermer and Ramaty 1986). 

We consider neutron and antineutron production In a relativistic plasma 
using a specific model with P co j = p co =500 MeV/c for all i. A low-momentum 
cutoff of 500 MeV/c corresponds to a low-energy cutoff of -125 MeV for the 
protons. We find that for these parameters and normalizations, the integrated 
production rates of neutrons and antineutrons are 3.2xl0~ 16 cm” 3 -s" 1 and 
5.5xl0 -19 cnf 3 -s _1 , respectively. This implies that -580 neutrons are 
produced for each anti neutron. 

OBSERVABLE CONSEQUENCES 

A review of the cosmic ray antiproton observations and models is given by 
Golden (1984), and the measured fluxes of antiprotons in the cosmic radiation 
are shown in Fig. 1. The prediction of the leaky-box model (Tan and Ng 
1983b), in which cosmic rays traverse an energy-dependent average grammage 
before escaping from the galaxy, is also shown here. This model fails by 
about three to five to account for the measurements between 2 and 10 GeV. It 
fails by over an order-of-magnitude to explain the low-energy measurement of 
Buffington et al. (1981), even if the effects of the solar modulation of 
cosmic-ray protons and antiprotons are taken into account (Tan and Ng 1983b). 

We propose that the antiprotons in the cosmic radiation in excess of the 
leaky box contribution have their origin in additional discrete galactic 
sources which produce antineutrons through nuclear reactions in an accretion 
plasma. The production of antiprotons through secondary collisions is 
suggested by the lack of observations of A > 1 anti nuclei, whose production in 
collisions of normal matter particles is negligible. This argues in favor of 
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Fig. 1. The observed 
cosmic-ray anti proton 
fluxes, shown by the data 
points, are designated by 
Bu (Buffington, et al. 
1981), Bo (Bogomolov, et 
al. 1979), and G (Golden et 
al., 1979). The lower 
curve gives the leaky-box 
model prediction of Tan and 
Ng (1983b) and the upper 
curve is the anti proton 
flux resulting from the 
decay of anti neutrons 
produced in a nonthermal 
relativistic plasma with 
p co = 500 MeV/c. This flux 
is fit to the anti proton 
data in excess of the 
leaky-box model 
contribution. 


an additional secondary source of antiprotons. Low-energy antiprotons are 
produced in relativistic plasmas as a direct consequence of the kinematics of 
such systems, and these plasmas could be formed in the vicinity of a compact 
object. The associated production of neutrons is a necessary consequence of 
this model, and suggests a potentially observable signature, which we discuss 
below. 


The antiproton flux resulting from the decay of antineutrons produced in 
a relativistic plasma with p co = 500 MeV/c and normalized to the data is shown 
in Fig. 1. The galactic antiproton flux was evaluated from the antineutron 
source spectra using the relation 
n (E ) m 


J (E ) = n 
P P 


_n n_ 

4irn„ 


(l 


. inelv' 
+ a _ ) 


( 8 ) 


esc 
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where ng [cnf 3 ] is the average proton density in the galaxy. The parameters 
describing the average grammage traversal A esc and the anti proton inelastic 
cross section o^ e1 are given by Tan and Ng (1983b). The effects of 
ionization energy loss and reinjection of antiprotons following non- 
annihilation inelastic interactions are unimportant for this source spectrum 
and are neglected. The constant n in equation (8) was chosen to give the 
proper source strength to account for the observed anti proton observations. 

We obtain a total galactic antineutron production rate of ~1.3xl0 39 V 67 n o< 2 
antineutrons per second for this model, where V 67 is the volume of the Galaxy 
in units of 10 6 cm 3 and n 0> 2 is the average galactic proton density in units 
of 0.2 protons per cm 3 . From Fig. 1 we see that the combined source of 
antiprotons resulting from this model and the leaky-box model gives an 
adequate representation of the data above 2 GeV. There is also a substantial 
flux of low energy antiprotons produced in this model, which might be 
sufficient to account for the low-energy measurement of Buffington et al. 
(1983), depending on the uncertain degree of solar modulation occurring during 
this period of observation. 

The value of n chosen to fit the cosmic-ray antiproton observations 
implies a related proton flux from neutron decay. We have evaluated this flux 
from the neutron source spectrum associated with the antineutron source 
spectrum, but now with oi in equation (8) replaced by the proton-proton 
inelastic cross section. The result is compared in Fig. 2 to the observed 
cosmic ray proton flux (data referenced in Gloeckler and Jokipii 1967 and Tan 
and Ng 1983b). As can be seen, the neutron-decay proton flux implied by this 
model does not conflict with the observed protons in the cosmic radiation. 
Above several GeV, these protons comprise about 10% of the observed cosmic-ray 
protons. At lower energies this source appears to contribute a larger 
percentage, but the effects of solar modulation will lower the calculated flux 
at these energies. The total galactic neutron production corresponding to the 
antineutron flux of Fig. 1 is -7.5x10 V 67 n 0 . 2 neutrons per second. 

Antineutrons and neutrons escaping from the interaction region either 
decay in interstellar space or are captured or annihilate in a dense medium, 
such as the atmosphere of a star. In the first case, positrons and electrons 
are formed in the decay of the anti neutrons and neutrons, but represent 
negligible additions to the electron and positron fluxes of the cosmic 
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Fig. 2. The data points 
give the observed solar 
minimum cosmic-ray proton 
fluxes (see Jokipii and 
Gloeckler 1967 and Tan and 
Ng 1983b). The solid curve 
represents the interstellar 
neutron-decay proton flux 
produced in association 
with the antiproton flux 
of Fig. 1. 


radiation due to the relatively small Lorentz factors and intensities of the 
neutron and antineutron source functions. An inner bremmstrahlung (Petrosian 
and Ramaty 1972) X-ray emissivity from the decay of the neutral particles is 
also produced at a level of - a-0.5 MeV/n • 10 42 n/s ~ 6xl0 33 ergs/s, where a is 
the fine structure constant. This represents a weak X-ray luminosity, even if 
concentrated in a single source. 

The capture of the neutrons and antineutrons is possible if, for example, 
the compact object is accompanied by a binary companion which subtends a 
significant solid angle about the interaction region. This companion star 
could also be the source of the accretion plasma that powers the system. 
Neutrons impinging on the atmosphere of the star would thermal ize and, 
depending on the 3 He content of the stellar atmosphere, produce deuterium in 
the capture reaction n + p -► 2 H + y(2.223 MeV), as in the case of solar flares 
(Wang and Ramaty 1974). If all neutrons are captured and produce deuterium. 
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this represents a production rate of ~3xl0 17 sec • 10 42 n/s ~ 3xl0 59 2 H produced 
during the age of the galaxy. This can be compared to a value of ~1.2xl0 63 c_ 5 
deuterons in the galaxy, where c_ 5 is the relative deuteri um-to-hydrogen 
galactic abundance in units of 1CT 5 . This source obviously does not conflict 
with the observed galactic 2 H abundance, but neither does it significantly 
contribute to it. 

The formation of a deuteron is accompanied by the production of a 2.223 
MeV gamma ray. For solar flares, where this effect is clearly seen, 
approximately 0.1 line photons are produced per neutron (Murphy and Ramaty 
1985). Let 0.1f_j be a numerical factor reflecting the probability that for 
each neutron formed in the accretion plasma, a 2.223 MeV photon is emitted 
into interstellar space. This factor incorporates a number of uncertain 
factors including the solid angle of the companion star and the probability 
that the gamma ray is both produced and able to leave the stellar atmosphere 
without being absorbed. We therefore have that ~ 7.5xl0 4 °f_iV 67 n o> 2 line 
photons per second are produced in the galaxy. If there are a large number of 
sources distributed in the galaxy, this represents a flux observed at Earth of 
- 10' 5 f_ 1 V 67 n 0<2 photons/(cm 2 -s-rad), normalizing to the observed flux of 
1.809 MeV photons from 26 A1 decay (i.e., 3 M 0 of 26 A1 with a mean life t 2 6 A 1 ~ 
l.lxlO 6 yr produce about 5xlCT 4 photons-cm' 2 -s _1 -rad -1 [Mahoney et al. 1984]). 

If the number of discrete sources is small, then the 2.223 MeV line flux 
could be observed from a single object. Taking the number of sources in the 
galaxy to be 10N lo , the flux of 2.223 MeV photons observed at Earth from a 
single source is ~ 6xlO -7 f_ 1 V 67 n 0# 2 /(N 10 r 10 2 ) photons/(cm 2 -s) , where ri 0 is 
the distance of the source to the Earth in units of 10 kpc. The detection of 
the line flux would be aided by the fact that the line intensity will vary 
with the underlying binary period of the system. In addition, since this 
estimate is based on the average source luminosity over the age of the cosmic 
rays, a temporary outburst could increase the emission expected from a single 
source. Likely candidates for observation include the binary X-ray sources 
listed in the Introduction, whose ultra-high-energy gamma-ray flux is 
evidently powered by an accretion mechanism, perhaps in association with a 
neutron star. A current interpretation of the origin of the high-energy 
gamma-ray flux pictures energetic hadrons emerging from the compact object and 
impinging upon the atmosphere of a companion star (e.g., Berezinsky 1977; 
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Vestrand and Eichler 1982; Hi 11 as 1984b; Kazanas and Ellison 1986). The 
observation of a 2.223 MeV signal in the transient event of 10 June 1974 
(Jacobson et al. 1978; Ling et al. 1982) is perhaps an example of this type of 
source, although with nonsteady accretion. 

CONCLUSIONS 

We have presented the physics of secondary neutron and antineutron 
production from particle collisions in a thick target and in a nonthermal 
relativistic plasma. Because the neutrons and antineutrons are not confined 
by the magnetic field of the accretion plasma, they can leave the interaction 
site and either decay in interstellar space or be captured on a companion 
star. Assuming that most of the neutral particles decay in interstellar 
space, we have fit the spectrum of secondary antineutrons to the observed 
cosmic ray antiproton flux, using a relativistic plasma with a low-momentum 
cutoff of 500 MeV/c and taking into account the leaky-box contribution of 
antiprotons. This normalization implies a production rate of about 10 42 
neutrons per second in the galaxy. This level of production does not violate 
the deuterium content of the galaxy and provides only a weak X-ray source 
resulting from inner bremsstrahlung in neutron decay. If about one 2.223 MeV 
gamma-ray photon is emitted into interstellar space for every 10 neutrons 
formed, for example by neutron capture on protons in the atmosphere of a 
stellar companion of the compact object, a flux of line photons of ~10 _s 
photons/(cm 2 -s-rad) could result from distributed sources in the galaxy. A 
gamma-ray line flux at 2.223 MeV might be observable from discrete sources 
depending on source luminosity, distance and duration of the active phase. 
Detection of this line could be facilitated by the variation of the flux with 
the period of the binary system. 
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